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Abstract 

Various models leading to predictions of negative capacitance, C, are 
briefly reviewed. Their relation to the nature of electric control is dis- 
cussed. We reconfirm that the calculated double layer capacitance can 
be negative under cr-control - an artificial construct that requires uniform 
distribution of the electrode surface charge density, a. For instance. It is 
shown that the combined relaxation of the ionic and electronic contribu- 
tions can result in C < even for the local statistical ionic models with 
strictly positive diffuse layer capacitance. In reality, however, only the 
total charge q (or the average surface charge density a) can be experi- 
mentally fixed in isolated cell studies ((/-control). For those a where C 
becomes negative under cr-control, the transition to g-control (i.e. relaxing 
the lateral change density distribution, fixing its mean value to a) leads to 
instability of the uniform distribution and a transition to a non-uniform 
phase. As an illustration, a "membrane capacitor" model is discussed. 
This exactly solvable model, allowing for both uniform and inhomoge- 
neous relaxation of the electrical double layer, helps to demonstrate both 
the onset and some important features of the instability. Possibilities for 
further development are discussed briefly. 



1 Introduction 



The admissible sign of the differential capacitance at charged interfaces was 
questioned in early 1970s in relation to the so called Cooper-Harrison catas- 
trophe |S], an apparent prediction of C < for "dipolar capacitors". Soon 
thereafter, similar anomalies were predicted for some ionic models [3 114) and, 
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somewhat later, for microscopic "relaxing gap capacitor" (RGC) models which 
accounted for the metal electron contributions in the interfacial capacitance 
(see E7| 1^ and references therein) . It is generally accepted that C 
must be strictly positive for open " (/i-controUed" systems, where the electric cell 
is connected to a source of controllable voltage ^Hl- However, it was shown 
that C < is possible as a stable state of an isolated RGC if electric control is 
maintained by a macroscopically uniformly distributed surface charge density 
a, so called cr-control (see ^11211221 for review). Prior to this discovery, an at- 
tempt had been made to prove, on general statistical-mechanical grounds, that 
for an equilibrium cr-controUed system C must be strictly positive 0. Some- 
what later improved analysis [2^1 UHl 1211 01 1^ showed the model ^ does not 
forbid negative C values. For both ionic and "relaxing gap capacitor" models, 
which together cover a very wide range of interfacial phenomena, the equilib- 
rium capacitance under cr-control can be negative. This is true for both the 
"compact layer" capacitance and the total double layer capacitance, including 
the "diffuse" layer contributions '24', "29'. 

Here we show, that though acceptable for cr-control, negative EDL capac- 
itance is not possible for a real isolated system, where uniformity of a is not 
enforceable and only the total electrode charge can be fixed. In other words, 
we show that in a domains where C is negative under cr-control, the system is 
unstable with respect to transition to an inhomogeneous state with nonuniform 
lateral distributions of the electrode charge density and of mobile electrolyte 
ions. To demonstrate this anomaly, we discuss a transparent "membrane ca- 
pacitor" model, which exhibits a C < domain presuming uniformity while 
becoming unstable in this domain if non-uniform surface charge distributions 
and membrane deformations are considered. 

This result addresses questions raised recently ^1 |S] with respect to the 
meaning and physical reality of C < for the isolated capacitor. Our analysis 
suggests that a model demonstrating this anomaly under cr-control can be used 
to study transition to an inhomogeneous interfacial state under g-control. 

2 Control of electrified interfaces - theory and 
experiment 

1. Experimental study of the electrical double layer (EDL) at electrochemical in- 
terfaces is usually conducted under "potential" {(j)-) control, where electrodes are 
connected to a voltage source. Changing the applied voltage in increments A(/) 
leads to corresponding changes of the electrode charge, Ag. Similarly, controlled 
modulation of the voltage, (p = 4>{'^)i results in charge modulation measurable 
by impedance techniques. Connection to a potentiostat, required to maintain 
(/i-control, results in an open system, which is treated by grand canonical meth- 
ods. 

The differential capacitance (per unit area) is then defined by the derivative 

C0 = d^a (1) 
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where the average surface charge density is a — q/A and A is the the surface 
area of the electrode. The form of Eq. ^is a typical response function 



(2) 



where F is the external parameter ("force") and X the conjugate intensive 
variable. 

2. The electrical properties of interfaces can equally well be studied by 
controlling the electrode charge, q. In g-control the measured quantity is the 
corresponding potential (f). Charge can be regulated by connecting the electrodes 
to a battery for brief periods of time At, measuring the current j, with the charge 
increments found by Aq = jAt. With q fixed, the resultant (j) is measured in an 
isolated system, i.e. a canonical ensemble. 

The corresponding response function, analogous to Eqs. ^ and |21 is the 
inverse differential capacitance 



Obviously, g-control is a synonym for a-control: fixing the total charge q is equiv- 
alent to fixing the average surface charge density a. The thermodynamic poten- 
tials of the expanded ((/)-controlled ) system, A^,a.nd the isolated (g-controUed) 
system, Ag^ are related by the Legendre transformation, 



3. Most calculations of the electric double layer (EDL) assume the electrode, 
often described as a charged fiat wall, has a uniform, fixed charge density. This 
is effectively tr-control, which is generally not equivalent to q- (or W-) control. 
The terms are interchangeable only if the equilibrium surface charge density 
is uniform on a scale exceeding atomic dimensions, i.e. if a = a = const. A 
counter-example is one where the optimized local charge density is non-uniform 
in the electrode plane, a = (j{rs) ^ a [ts is the radius vector in the electrode 
plane), reminiscent of "charge density wave" states in electron plasma. In prac- 
tice there is no way to control the local charge density; in isolated systems only 
the total charge can be constrained externally, tr-control is a purely theoretical 
construct; its predictions must be tested to determine if a uniform surface charge 
density and its corollary, a laterally-uniform ionic distribution correspond to a 
real equilibrium state. 

3 Admissible sign of the differential capacitance 

3.1 (/)-control: the open system 

General thermodynamic J7| and statistical-mechanical JHI treatments of elec- 
trified interfaces show that, under 0-control, differential capacitance must be 
strictly positive. In our view attempts to circumvent this restriction pi I35L 
have been based on misinterpretations of the nature of potential control j27ll29| . 



(3) 



A^{(J,(t)) = Aq[a) ~ q<j). 



(4) 



3 



Near a critical voltage </)cr, defined by C~^{(j)cr) — 0, the system becomes un- 
stable. The transition to a new state is accompanied by charge flow from the 
potentiostat to the electrodes, a sort of electrical " breakdown" (see ^fl ^1 1^ 
for more details), which would be a unique path to phase transformation as- 
suming lateral uniformity. However, as discussed below, the transition can also 
involve formation of a laterally non-uniform phase accompanied by nonuniform 
redistribution of the electrode charge density cr(p). 

As the requirement that > is now generally accepted, we turn to 
treating isolated systems. The sense of the upcoming discussion is already 
implicit, once having recognized that the first of these instabilities is forbidden 
by the very definition of "g-control" since, once q is fixed, electric contact with 
the potentiostat must be interrupted. Can C be negative under this constraint? 
As the admissible sign of C in isolated systems has almost always been analyzed 
in (T-control terms, we first consider this case and postpone discussion of the 
more general g-control. 



3.2 cr-control in the isolated system 
3.2.1 Primitive models of electrolytes 

Interest in the admissible sign of C in the theory of the diffusive layer was 
stimulated by work of Blum, Lcbowitz and Henderson j4j . They tried to provide 
a rigorous restriction on the sign of C for "primitive ionic models:" charged hard 
ions in a uniform dielectric medium between two rigid, uniformly charged walls. 
The corresponding Hamiltonian is quite generally |26| 

{R}) = ^ - ^/({R}) + H ({R}) (5) 

where {R} refers to a particular configuration of the system (charge coordi- 
nates, dipolar orientations, etc.). The first term describes direct interaction 
between the charged walls, with d the inter-wall distance. (t/({R}) accounts 
for interaction between the electrolyte and the electrode field (the physical sig- 
nificance of / will be clear shortly) and H' is a a independent interaction energy. 
The potential drop between the charged plates is 

^(a) + <f> (6) 

eeo 

where 

_ /e^^^("'W)(...)rfO 

is a canonical average with integration over the system's configurational space, 
n, and /3 — 1/kT. Eq. reveals the meaning of /: < / > is the potential 
drop induced in the electrolyte by the field of the charged plates. It arises 
from redistribution of free (ionic) charges shielding the applied field, and from 
repositioning of the bound charges ( the reorientation of molecular dipoles). 
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For Hamiltonians of the type of Eq. |31the capacitance satisfies the general 
condition j2B| 

leading to a simple and self-evident result: 

< — . (8) 

eeo 

Eq. |S1 simply means that redistribution of free charges and molecular polar- 
ization induced by the electric field in the electrolyte reduce the potential drop 
between the electrodes and increase the capacitance C, a condition that places 
no constraints on the sign of C under cr-control. One should note here that if 
the distance d >> Au, the characteristic Debye length in the electrolyte, then 
the total inverse capacitance of the cell splits into two independent double layer 
contributions belonging to two "electrodes:" 

If the sign of the total cell capacitance C is unrestricted under a-control, this is 
even more true for the individual double layer contributions, Ci and C2. 

Eq. [7| (Eq. 24 of (2^1) was derived in a study of a "dipolar capacitor" 
("DC"), a lattice of point dipoles embedded between the plates of a parallel- 
plate capacitor, a model often used for analyzing a compact layer at metal- 
solvent interfaces. For the DC 

is the potential drop corresponding to an arbitrary configuration of the molec- 
ular dipoles with Pz the average surface density of the dipole moment in the 
lattice and pi^z the projection of the individual dipole moment normal to the 
surface of the lattice. 

If interaction between the charged walls (the first contribution to the Hamil- 
tonian, Eq. [S]) were properly included in the Hamiltonian of the primitive ionic 
model 0], it would also lead to Eq. [7|(see p. 68 of 1201) with 



/ — fion / 



qiZi 

eo ^ 

where qi is the charge of i-th ion and Zi its distance from the charged wall po- 
sitioned at 2; = 0. Eq. [71 and its analogs have been repeatedly derived and 
discussed j^SlElEIl; they hold for any model in which the electrodes are treated 
as hard charged walls with distinctly separate regions occupied by electrode and 
electrolyte. With these restrictions the interaction between the electrode and 
the electrolyte can quite generally be described by a contribution ~ (J^qiZi 
where the summation includes both the mobile ionic and the molecular multi- 
pole charges j^Zl- Such constraints clearly exclude " polarizable" models, those 
explicitly treating molecular electronic polarizability, electron density penetra- 
tion into regions occupied by electrolyte, etc. 
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3.2.2 Relaxing gap capacitors (RGC) 

Immobility of the charged "plates" in primitive models does not account for an- 
other important phenomenon, possible displacement of the "electronic plate" of 
intcrfacial capacitors and of the equilibrium positions of the electrolyte species 
in contact with the electrode, in response to charging |271 1^ EHI ■ These ef- 
fects are effectively illustrated by the "relaxing gap capacitor" metaphor, which 
emphasizes the dependence of the effective gap d on charging. The potential 
drop can be quite generally represented as 

(j) = —ad{a) (9) 

££0 

where d{a) is the effective separation between the "plates" of the capacitor 
associated with the "centers of mass" of two microscopic charge distributions 
(see 1^ ^1 1^ for review) . We assume a uniform dielectric background e 
between the plates, in the volume occupied by free charges. For models using a 
non-uniform background, typical of unified models used to account for both the 
" inner" (or Helmholtz) and the diffuse layers , the expression is more complex 
and must also account for the distributions of the bound (polarization) charges. 
However, these details are not essential; they are omitted here. 
For the electrode-electrolyte interface the effective separation is 

d{a) ^ Zi{a) ~ Zi^e{<^) (10) 

where 

_ / Plezdz 
a 

with z the coordinate normal to the electrode surface. The inverse differential 
capacitance for the RGC is 

( — C)-i ^ iad{a)y, = d{a) + ad(a);. (11) 

eeo 

Dependence of d on ct, a general feature of practically all double layer models, 
imphes that C is dependent on a (or on the apphed voltage). The "plate" 
displacement that contributes to the variation of d not only reflects a shift of 
charge density profiles, but is more a consequence of shape variation |29j . Quite 
typically, there is always a a range in which charging decreases the effective 
gap. Elastic compression of the lipid membrane by electric stress |51ES1|2H1 and 
response of the Gouy-Chapman-Stern (GCS) diffuse layer to charging are 
two representative examples. 

In a range of a where the effective gap contracts with charging, d (a) < 
and Ca can be negative if 

Aa -d- (J ■ Ad<0. (12) 

This inequality means that a potential increase due to a change of a is over- 
whelmed by its decrease due to gap contraction. A number of electromechanical 
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and microscopic models |^1|27||^|2HI show that negative capacitance (NC) un- 
der (T-control is compatible with system stability. For instance, for every fixed a 
(including the domain where C < 0) the equilibrium gap of the elastic capacitor 
is defined by a stable balance between elastic and electrostatic forces. Similarly, 
density functional minimization led to a NC at metal-electrolyte interfaces due 
to relaxation of the "electronic plate" of the capacitor [111 1121 [TC] . 

Previous work |12[ 1^ showed that accounting for " electronic plate" relax- 
ation in combination with traditional GCS and similar models typically leads 
to negative C domains even though the GCS model itself (as with any other 
"local" statistical model where ionic concentrations are local functions of the 
potential) does not lead to such anomalies a point further demonstrated 
in the Appendix. Thus, while we agree that finding C < under u-control 
must be both common and important for ionic models of electrolytes |13| . it is 
not a necessary condition for the appearance of this anomaly. If other relax- 
ation mechanisms are taken in account a NC domain may arise even if the ionic 
contribution is positive. 

3.2.3 Possible capacitance anomalies for local u-controlled statistical 
models with interfacial relaxation 

Consider a conventional two-layer model of the double layer, with a Helmholtz 
layer accounting for the finite ion-electrode distance of closest approach, a, and 
the diffuse layer accounting for the electrolyte's ionic charge distribution [HlEl. 
Its inverse capacitance is 

C-\a)^CH\a)+Co\a) (13) 

Here 

CH\<T)^d^^H (14) 

with (f>H and the potential drops in the Helmholtz and diffuse layers respec- 
tively. We now treat "local statistical" diffuse layer models. 

These simple models describe the ionic density pi at a distance z from the 
electrode as a function of the local potential (p{z): 

Pi = PiMz)] (15) 

The classic example is the Poisson-Boltzman-Gouy-Chapman (PBGC) model of 
symmetric electrolytes 

Pi — qriQ cos h[~l3q(p{z)] (16) 

with no the bulk concentration of cations or anions. For the slightly more com- 
plex model treating "ionic saturation," finite ion size driven entropic restriction 
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on local ionic concentration, which can be especially significant in solid elec- 
trolytes, the ionic density is 

f^_ f exp[-/jy)] exp[(3ip] 

4-ei + eiexp[-/3^] i-e2 + e2exp[-/3(^] 

where 0i_2 = ?^o/-^i,2 and are maximal possible concentrations for the 
cations and anions respectively. Concentration limitations reflect finite ionic 
size and, for solid electrolytes, the limitations on the number of possible ionic 
defect sites in the crystal lattice. 

The one-dimensional contact model is illustrated in Fig. 1. 



Fig. 1 Schematic representation of the ionic, pf, and the electronic, 
/9g , charge distributions induced by the electrode surface charge 
density a. 

We use a two-layer dielectric model of the interface 

where Sh and are dielectric constants of the Helmholtz layer and bulk 
electrolyte respectively. Ignoring the penetration of the electron density into 
the "diffuse" layer the potential drop in the Helmholtz layer is 

= -^dnia) (19) 

where 

dH{<7)=a-Zeia) (20) 
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is the effective gap of tlie Helmlioltz layer. As a result, in direct analogy with 
Eq.ini 

Ch^{<^) = [a - -a da Ze] (21) 

£h£q 

The significance of Eq. Elis that dif (cr) depends on a via displacement of the 
"electronic plate" Ze of the equivalent interfacial capacitor during the charging 
process. This effectively turns Ch into a "relaxing gap capacitor" [5^11^1^ . 
The cr-dependence of the position of the equivalent electronic plate has been 
studied using a density functional approach for various models of the ionic 
charge distribution near the electrode, from the Raleigh picture where the 
ionic charge density is localized in one or two monolayers nearest the electrode, 
to the "external field" limit, where countercharges are located far from the 
electrode )■ The diffuse layer distribution is generally bounded by 

these limits; Ze{<j) is always well approximated by the cubic polynomial 

Zg(cr) = Ze{0) + s a + pa^ + ra^ (22) 



= ^—{a- Ze{0) -2s cr + 3pa^) (23) 
The diffuse layer ionic distribution is related to a via the strict "sum rule" 

= -2e„£o / p{'p)d(p (24) 
Jo 

which leads to a general expression for the capacitance 



£v£o Pi{4>D) 

where pi{(f)D) — Pi{a) (see Eq. 1151 ) is the ionic charge density at z = a, the 
point where the local potential tp equals (po (we choose <p(cxd) = 0). Since a and 
Pi (a) are of opposite sign, this is simply 



Pi[<i>D) 



(25) 



which shows that Cd is always positive and finite for any finite a. Thus local 
statistical models don't satisfy the criteria suggested by ^31 ■ We now show that 
charging induced relaxation of the "electron plate" (see Eqs. 1221 and VF^ can 
lead to a negative total capacitance even if the diffusive contribution is strictly 
positive. In the spirit of a-control we express (po and Go solely in terms of a. 
Corresponding relations were discussed for various local ionic models in |3(J[I15| . 
Thus for the model of Eq. El the result is 

(1 - Gi + eiu)(l - Oi + QivT^Y = G{a) (26) 
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where 



/ -i-PD) \ r<f \ I ^ 2\ Qi N2 

— Vl 



/cT / ' ' ' ^ \2eeoNikT J 62 iV] 



For the case Ni = N2 = N, corresponding to a solvent with cations and 
anions of equal solvation shell radii or solid electrolytes with Shottky defects, 
this can be solved analytically. The capacitance is 

Cd =eeono . . + 4^(1 - 60) (27) 

where S ^ {G - IJ/Qq, Qq = uq/N. For A^i ^ A^2, Eq. Eimust be solved 
numerically. We limit consideration to comparatively low ionic concentrations 
where restrictions on ionic packing can be neglected and the ionic distribution 
is described by the PBGC model, Eq. ^1 Then Eq. [pleads to the familiar 
expression 

cr^^ = «(/?+ ^Y^' = + i'y/'^on/' (28) 

with 



Ld = xl-TTT^ = 2-82 IQ-^^x —[m], ao = ^J^n^esokT = 7.67 10" VcqeT [C/ 
V 2(7-^no V Co 

(29) 

where Cq is the molarity of the solvent. In the a — Q limit the capacitance is 
naturally determined by the electrolyte's Debye length. 

To explore the analogy with the relaxing gap capacitor, it is instructive to 
represent (C|^^'^'^)^^as ~dcr{(Tl{<y) ) and determine the cr— dependence of the 
effective gap I. Solving 



Ld 



[l + (a/ao)2]i/2 
yields 

1(5)^^ Log E + + (30) 

where S — a/ao. The dimensionless gap relaxation a = I jLjj (line 1) and the 
corresponding dimensionless potential Eq;(S) (line 2) are shown in Fig. 2 and 
compared with similar features of the elastic capacitor ( lines 3 and 4). 



Fig. 2. Charge dependence of the dimensionless effective gap a(S) 
(curves 1 and 3) and the potential $ — Sa(E) (curves 2 and 4) 
for the PBGC diffuse layer (curves 1 and 2) and the membrane 
capacitor (curves 3 and 4); S is the dimensionless charge density. 
For the PBGC diffuse layer a{T.) = 1{T.)/Ld (see Eq. OOland the 
definition below) and for the membrane capacitor a(E) = 1 — 
(compare with Eq. 1331 below) 
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s (dimensionless) 
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This picture clearly demonstrates how the relaxing gap capacitor is related 
to the GC diffuse layer. In both models charging contracts the effective gap. 
However, the rate of contraction in the PBGC model is insufhcient for the 
formation of a peak in (curve 2), observed in the elastic capacitor (curve 
4) and associated with the C < domain. This follows directly from a general 
property of all local models, that C is positive at all surface charge densities, a. 

By considering the solid electrolyte, AgCl, for which the compact layer elec- 
tronic properties have been studied in ^1 130] , we show that a compact layer 
can lead to a negative total (diffuse + compact) layer capacitance. The tem- 
perature dependence of the parameters e and cq are given by PHI- At 400°C, 
Co ~ 43 mM and e 3 so that Ld 2.7 A and cto ~ 0.55 /iC/cm^. Combining 
Eqs. ESlandEHlwefind: 

^0^"' = -JT-r^r^anFi + — (a - a;e(0) - 2s a - ipa^ - ^rs^) (31) 

According to the commonly held view, the contribution of the compact layer 
effectively increases the effective gap, thus increasing the total inverse capaci- 
tance. However, the appearance of a domain where C < implies that C~^can 
be reduced sufficiently to become negative. Therefore , one expects that if local 
model diffuse layer contribution is strictly positive, this should be even more 
pronounced in presence of the compact layer. The following discussion demon- 
strate this to be wrong. Results for the double layer capacitance are presented 
in Fig. 3. 

The traditional model (line 1) includes the diffuse layer, Eq. |2H1 in series 
with the fixed compact layer contribution calculated in the " perfect conductor" 
model. This ignores diffuseness of the electronic distribution and its relaxation 
in the compact layer, in effect assuming 

^-1 ^ £o£g 
" a 

Naturally this capacitance is positive at all a. Accounting for electronic re- 
laxation in the compact layer dramatically affects behavior. The parameters 
Xe{o), s, p and r (Eqs and |23J| were calculated by ^l|3n| for Au/AgX con- 
tacts. 

Fig. 3. Double layer capacitance C for two models of the interface:l 
- traditional PBGC-Helmholtz model (the position of the "effective 
electronic plate" is fixed at = 0); 2 - modified PBGC-Helmholtz 
model, which also accounts for relaxation of the effective "electronic 
plate." 



The corresponding C is shown for AgCl (line 2). Electronic relaxation leads 
to a pronounced capacitance asymmetry and steep growth in the cathodic range 
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of charges, both eflFects observed experimentally [221 E3|- It yields a vertical C 
asymptote separating C > (a > acr) and C < (not shown) domains. This 
result is general, and unaffected by variation of er (typically from 2 to 5), e„ 
(from 3 to 10) or electronic parameters reflecting different approximations to 
the electronic density functional ^1 OOl E]- The appearance of a negative 
capacitance domain is similar to phenomena arising from the compact layer at 
metal-solvent interfaces |16| Thus, C < may occur under (t— control even 
for local diffuse layer models if the metal electrode modeled realistically. Put 
differently, a C < domain may arise even for local ionic models with strictly 
positive diffuse layer contributions, if additional electronic (and possibly other) 
mechanisms for interfacial relaxation are considered. 



3.3 q-control in the isolated system 

As discussed, cr-control is an artificial construct. Nevertheless, charges are nor- 
mally distributed uniformly in the plane of the electrode and the charge q nat- 
urally yields a uniform charge density cr, in which case a- and g-control are 
identical. However, it is possible that under special conditions a non-uniform 
distribution of charge in the electrode plane becomes energetically preferable. 
Under these conditions the cr-control metaphor implies that artificial restraints 
(forces) were applied to the surface charges to enforce uniformity of a. Elimina- 
tion of these artificial restraints would result in transition to an inhomogeneous 
state. We will show this and its relation to the sign of the capacitance by 
considering an undulating membrane capacitor, an exactly solvable model that 
contains many features common to real EDLs, which is thus generally useful for 
discussing double layer behavior under critical conditions. 

For a membrane capacitor under potential control the onset of instability 
results from the steep increase of the electrostatic force ~ (fp' /h? where (f) is 
the applied potential and h is the membrane thickness. This arises because 
the system is open: as h decreases, thinning leads to charge transfer between a 
battery and the plates of the capacitor required to maintain the fixed value of the 
potential drop (f> ^ ah — const. Mathematically, charge transfer is controlled by 
the term — relating the thermodynamic potentials of the open and isolated 
capacitor (Eq. ^J. Thus, both the charge density and the attractive force 
increase unrestrainedly as the plates approach one another. At a critical voltage, 
no increase of the elastic repulsive force can compensate for the catastrophic 
increase of the electrostatic attraction, leading to the potential-driven instability 
first noted by Crowley j^. 

In an isolated system this pathway to instability is forbidden. When isolated, 
the total plate charge is fixed and uniform system thinning cannot increase the 
attractive force, which remains constant, ^ cr^. Instability can only result from 
nonuniform redistribution of the charge density in the plane of the membrane 
capacitor with an associated inhomogeneity in h. Thus we are led to consider 
the possible lateral instability of an isolated flexible membrane capacitor. For 
</)-control we solved this problem for electrolyte charge distributions that satisfy 
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the Poisson-Boltzmann equation Here we consider a simplified example, a 
capacitor in contact with a "perfect conductor," i.e. Ssoivent — oo, so that the 
membrane surfaces are isopotentials, a reliable approximation even for dilute 
electrolytes pS] . 

Consider the parallel-plate membrane capacitor with its midplane at z = 
and an unperturbed thickness (when q — Q) equal to /iq. Charging the membrane 
to charge density a yields an electrostatic pressure which gives rise to membrane 
compression. Introducing the thinning coefficient a — h/hg, with h the thickness 
of the compressed membrane, the total energy of the uniform slab is 

Wo = Wdfi + We,o where Wd,o = lrKs{a - if and We,o = x^^^^; (32) 

2 Zeeo 

the two terms are the harmonic approximation to the deformation (stretching- 
compression) energy and the electrostatic energy respectively. Kg is the stretch- 
ing modulus and the index "0" refers to a uniformly deformed membrane. The 
equilibrium membrane thickness h{a) — hoa{a) is found from the condition 

daW = 

leading to the thinning coefficient 

a{a) = l-h — f (33) 

O (Tcr 

where 



which is interpreted in what follows. The transmembrane potential drop and 
the corresponding inverse differential capacitance are 

4>= ah{c7) = —hoa[l ~ ^{a/acrY] (35) 

£So SSq o 

and 

C-\a)^^^C,'[l~ {a/a^rn (36) 



where 

''-■0 



Cn ^ — ho 



eso 

is the inverse capacitance of a capacitor with the fixed gap ho- Eq. 1361 reveals 
the meaning of Ucr- It is the charge density where the differential capacitance 
Cfj becomes infinite; Co- is negative for \a\ > a^r- The corresponding membrane 
thickness is 

2 

at (7 = (Tcr the membrane has thinned by 33% , a value typical of " relaxing 
gap" capacitor models |31l 1161 l?7l 1241 OKI I29j . The membrane is stable relative to 
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virtual uniform compression (thinning) under cr-control (i.e. assuming a uniform 
surface charge density) for all a including the range where C < 0. This is 
verified from Eq. 

dlaW\a = — > 0. (37) 

a 

We now focus on the energy change, ISW , for an isolated membrane capacitor 
in response to a small charge increment, Act, assuming a-control. Using the 
relation 

we find: 

ISW = 0(ct)Act + ^C'-I(ct) (Act)2. (38) 

This equation has important consequences. First virtually separate the mem- 
brane capacitor into equal patches J and //, each of area Ajl and permit the 
charge Ag — A/S.a/2 to flow from / to //. For simplicity neglect boundary ef- 
fects and assume each charge density is uniform (cti = ct — Act and CT2 = ct-|-Act) 
and that the patches deform independently, i.e. two membrane capacitors are 
in parallel, and not elastically coupled. Since the potential is constant in the 
plane of the membrane, we find from Eq. 1381 that the total energy change is 

AM/12 = ^ AM^i + ^ AW2 = ^ C'-^ (Act)^ (39) 

Thus ATyi2 is negative if C < for the chosen ct. Put differently, C < pro- 
vides a driving force leading to a nonuniform charge distribution and membrane 
deformation if the artificial CT-control restrictions are relaxed. 

The energy penalty results from the continuous transition between the prop- 
erties of the two membrane patches and can be described as a linear interfacial 
tension. It is proportional to the length of the border between them and, for 
large A, is negligible relative to AW^i2, Eq. 133 Thus this charge density re- 
distribution and the corresponding non-uniform deformation of the membrane 
is possible energetically for those charge densities that lead to a C < domain 
assuming CT-control. The appearance of C < in treatments that presume a 
uniform charge density indicates the system is unstable . This result is similar 
to the thermodynamic arguments of Nikitas |21| who considered equilibrium 
conditions between two separate surface phases. We will now show that the 
prediction of negative capacitance under CT-control also implies that there is the 
possibility of forming an inhomogeneous phase under g-control.. In our virtual 
experiment we assumed the charge density is uniform in each patch. Releasing 
this restriction provides other pathways for transition to a non-uniform state. 
Consider, for example, membrane stability relative to symmetric undulations, 
the harmonic variation of membrane thickness: 

h{x) = h + 2 u cos{kx), 
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where u is the amplitude of the undulation of the membrane surfaces; the cor- 
responding "left" and "right" interfaces are described by the equations 

Zr,i{x) ~ ± zq{x) where zo{x) = h/2 +u cos^kx). 



This problem has been discussed previously for (/)-controlled systems (see |H1 
I25L I29j and references therein). We treat a ^-control environment by fixing the 
average charge density q — a A instead of the membrane potential (j). 
Some aspects of the solution procedure should be stressed. 

(1) Unlike under 0-control, under g-control the transmembrane potential 
drop V is not fixed by the external source (battery, potentiostat). However, the 
conductive surfaces are still equipotentials and the potential 0^ is constant on 
the membrane plane. 

(2) The value of (j)q depends on both the original charge density a of the 
unperturbed membrane and the parameters u and fc, characterizing the undu- 
lations. 

(3) (j)q is then determined as follows: 

(a) The solution for a fixed but arbitrary (j) determines the poten- 
tial v{x, z) within the membrane \25\ (see also discussion in ^ and references 
therein). 

(b) The equation A'Ka{x) = eV„w(x, zo(a;)) determines the local charge 
density cr(a;), where V„ is the normal derivative taken at the interface zq{x). 

(c) The total interfacial charge q{V), is found by integrating a{x) over 



the interface with a weighting factor, 1/^/1 — [dxZo{x)]'^, that accounts for the 
membrane stretching associated with undulations. 

(d) (pq is determined from the condition q{V) — a A: 

a , , fccoth(fc/i/2) n, 

0, H 1 ) (40) 

with h = h{a) = hoa{(7). The externally fixed potential in the equations 
for (/)-control can now be replaced by (jjq, which completes the solution of the 
problem for (/-control. 

The membrane's electrostatic energy is then 

= Wl, + Wl^ (41) 

where 

Wl,^^h{a) (42) 
is the energy of the uniform membrane slab and 

^e% = -l^^,o^coth[^] (43) 
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is the undulatory contribution. The onset of instabihty is determined by com- 
petition between the decrease of the electrostatic energy, Eq. 021 and the cor- 
responding increase in membrane deformation energy averaged in XY-plane, 
Wd,u- For simplicity, we consider small k (the long- wavelength limit) kh << 1. 
Similar to 23 where a slightly different form of Wd,o, Eq. 1321 was used, we can 
represent Wd,u as 



Wd,u + 0{{khf)]. (44) 

"-0 

Higher order terms in (fc/i)^ arise from surface tension and bending contributions 
to the elastic energy and are neglected. They are analogs to the non-uniform 
interfacial contributions of the previous example. In the same limit Eq. 1431 can 
be represented as 

Wl^ - -2Wiy[j^ + Oiikhf)]. (45) 

The uniform distribution becomes unstable when 

Wu = Wd.u + wi^ < 0. 
Substituting Eqs. 1431 and OSl we find the condition of instability: 

1^1 >0-cr. 

which is equivalent to 

c-i < 0. 

Our thought experiment presumed that the conditions for cr-control , i.e. unifor- 
mity of cr, could be arbitrarily relaxed anywhere within the negative C domain. 
In reality there is no way to enforce uniformity when the system is unstable; 
thus the distribution spontaneously becomes inhomogeneous at the edge of this 
domain with critical point a — Ucn where 

C-i(a,,)=0; (46) 

the transition actually occurs a bit earlier, at a point roughly determined by a 
Maxwell construction '27". 

We have analyzed a simplified model corresponding to a concentrated elec- 
trolyte (Debye length Ad — > 0). Further analysis based on our previous work 
shows that Eq. 0^1 also determines the onset of instability for finite |25| . 
This result is valuable on its own. As already indicated the Poisson-Boltzmann 
approximation and other local statistical models do not predict NC (see j^l 177] 
and references therein). Eqs. 1111 and IT^ show that for these models the rate 
of gap contraction with charging, d'{a), is always less than 1/|(t| and thus Ca 
is always positive. Consequently, this class of models would not satisfy the cri- 
teria suggested in linking model quality to the appearance of a Ccr < 
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domain. Our results indicate that adding another relaxation mechanism im- 
mediately leads to the appearance of a < domain and instability. This 
and our earlier discussion of the electronic models illustrates that anomalies are 
more typical than expected based on purely ionic models with immobile charged 
plates. 

In previous analysis |24| we considered two elastically coupled membrane 
capacitors, with the extra term in the deformational energy ~ a (ui — 712)^ 
accounting for the non-uniformity penalty (differential thinning of the patches). 
Depending on the coupling constant a, this system could exhibit a C < 
domain before transition to a nonuniform state. Our present discussion implies 
such a picture is unrealistic. In terms of the first model, the constant a must 
be proportional to the width of the transition region relative to the area of 
the patches A and thus can become infinitesimal if A is sufficiently large. In 
addition, as shown in the second example, the non-uniformity contributions 
become insignificant at small k which provides a reasonable pathway for the 
onset of instability. 

Finally, we reiterate the major difference between q- and 0-control for the on- 
set of instability. Under 0-control (an open system connected to a potentiostat) 
stability is lost simultaneously for both uniform deformation and undulations 
|51I2S]. In contrast, a g-controUed (isolated) system is always stable with respect 
to uniform deformation (see Eq. |^and the corresponding discussion) and only 
loses stability in transiting to a nonuniform state. An important consequence is 
that in a (/)-driven transition the original and final phases correspond to differ- 
ent values of q while under g-control only the local charge density can change 
provided its average value is fixed. 

4 Perspectives for further study 

It is by now well established that the capacitance can be negative for uniformly 
charged surfaces under the artificial conditions of cr-control. Many statistical 
ionic models have demonstrated this anomaly (see [HH El ^] and references 
therein). While in the 1980s many would have viewed such predictions as model 
faults, the most recent view represents a dramatic change of mind, and it 

is even suggested that the prediction of C < must be considered as a criterion 
validating an ionic model J3] rather than an imperfection. We do not disagree. 
In fact, it accords with our observation |^E1|271|^ that such predictions are 
typical of models jointly accounting for various contributions (electronic, ionic, 
etc.) to the charging induced relaxation of the effective gap of the interfacial 
capacitor. 

What is the physical significance of such predictions? Our analysis indicates 
that they imply an instability with respect to a transition to an inhomogeneous 
state. To analyze the consequences, the artificial assumption of a uniform sur- 
face charge distribution must be dropped, given that in a real metallic electrodes 
the electrons are free to move and thus surface charge density can become, at 
least in principle, laterally non-uniform. This non-uniformity can be even more 
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pronounced in soft media such as lipid bilayers. The value of the models predict- 
ing this anomaly, emphasized in |13| . is even greater since they are candidates 
for analyzing such instability. 

We must re-emphasize that in our usage the terms "instability" and "tran- 
sition" are not related to real interfacial critical phenomena. While the phase 
transition actually occurs before becomes negative, we permit the system 
to enter this domain by artificially maintaining cr-control. Relaxing the uni- 
formity constraint at any a within the C < domain leads immediately to 
a transition to an inhomogeneous state of fixed q, q — a A. Although the ini- 
tial (cr-controlled) state is artificial, the final stable inhomogeneous phase (if 
it exists) is real since the equilibrium state is unique. Thus, our approach is 
useful for testing and developing statistical models that describe both uniform 
(C > 0) and inhomogeneous (regions with C < under cr-control) phases. 

Predicting instability does not necessarily imply that the model describes 
the formation of a new stable inhomogeneous phase. As the transition can 
be accompanied by a substantial local increase of the charge density and a 
corresponding local increase of the ionic density in the EDL, a model must 
have a stabilization mechanism that interrupts the propagation of instability. 
This would permit formation of new stable phase before the condition of " ideal 
polarizability" is broken and interfacial charge transfer occurs, fonic size and 
correlation effects in the electrolyte must be important here. 

Although the condition C < 0, obtained for a primitive ionic model of the 
EDL, can be an important factor leading to instability, it is not a necessary 
condition. Even with ionic models that by themselves do not lead to instability 
(such as GCS model) the addition of other mechanisms of relaxation, such as 
a displacement of the " electronic plate" of the interfacial capacitor, can lead to 
Ccr < and thus trigger the instability |12[ . In other words, this anomaly 
should be even more commonplace than is implied by ionic model studies. 

Another important question is a comparison of phase transitions predicted 
for the open and the isolated systems. While in the first case the transition can 
be accompanied by charging the electrodes, in the second case the lateral varia- 
tion of charge keeps the total charge fixed. Finally, the inherent inhomogeneity 
(roughness) of an electrode surface (especially with respect to solid electrodes) 
must be considered. The influence of the roughness on the equilibrium proper- 
ties of EDL is well established > ^^'^ its possible effect on the surface phase 
transition can also be a promising field for further research. 
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